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ABSTRACT 



Context. Laboratory spectroscopic research plays a key role in the identification and analysis of interstellar ices and their structure. 
To date, a number of molecules have been positively identified in interstellar ices, either as pure, mixed or layered ice structures. 
Aims. Previous laboratory studies on H 2 0:CO ices have employed a 'mix and match' principle and describe qualitatively how ab- 
sorption bands behave for different physical conditions. The aim of this study is to quantitatively characterize the absorption bands of 
solid CO and H 2 0, both pure and in their binary mixtures, as a function of partner concentration and temperature. 
Methods. Laboratory measurements based on Fourier transform infrared transmission spectroscopy are performed on binary mixtures 
of H2O and CO ranging from 1:4 to 4:1. 

Results. A quantitative analysis of the band profiles and band strengths of H2O in CO ice, and vice versa, is presented and interpreted 
in terms of two models. The results show that a mutual interaction takes place between the two species in the solid, which alters the 
band positions and band strengths. It is found that the band strengths of the H 2 bulk stretch, bending and libration vibrational bands 
decrease linearly by a factor of up to 2 when the CO concentration is increased from to 80%. By contrast, the band strength of the 
free OH stretch increases linearly. The results are compared to a recently performed quantitative study on ILO:C02 ice mixtures. It 
is shown that for mixing ratios of 1 :0.5 ILO:X and higher, the H2O bending mode offers a good tracer to distinguish between CO2 or 
CO in H2O ice. Additionally, it is found that the band strength of the CO fundamental remains constant when the water concentration 
is increased in the ice. The integrated absorbance of the 2152 cm -1 CO feature, with respect to the total integrated CO absorption 
feature, is found to be a good indicator of the degree of mixing of CO in the H 2 0:CO laboratory ice system. From the change in the 
H 2 absorption band strength in laboratory ices upon mixing we conclude that astronomical water ice column densities on various 
lines of sight can be underestimated by up to 25% if significant amounts of CO and CO2 are mixed in. 

Key words. Astrochemistry, Line: profiles, Molecular data, Molecular processes, Methods: laboratory infrared, ISM: molecules, 
Infrared: Band strength, ISM, Spectroscopy: Solid State 



1. Introduction 

Water and carbon monoxide are common constituents in vast 
regions of space, both in the gas phase and in ices. Interstellar 
water ice was first identified in 1973 via a strong band at 3.05 
yt/m and unambiguously assig ned to water ice following compre- 
hensive laboratory work dMerrill et alJI 1976t iLeger et ai1I l979; 
Hag en et all 1 19791) . Meanwhile, it has become clear that H2O 
ice is the most abundant ice in space. The OH stretching mode 
at 3.05 jum and the H2O ben ding mode at 6.0 nm are detected 
in many lines of sight (e.g. 
1990; iMurakawa et alj [2000; 



Willne r et all 



Boogert et al 



200UlGibbetal . 2004; K nez et 1 



19821; iTanakaetal 
2000; iKeane et al 



OOgl 

al. 2005D and in many different 
environments, ranging from quiescent dark clouds to dense star 
forming regions a nd protoplanetary disks dWhittet et al.l fl988; 
iTanaka et alJH 994). It has been a long-standing problem that the 
intensity ratio of these two water bands in astrophysical obser- 
vations is substantially different from values derived from labo- 
ratory spectra of pure H2O ice. In recent years it has been pro- 
posed that this discrepancy may be due to contributions of other 
species, in particular more complex organic ices, to the overall 
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intensity of the 6 fim band dGibb & Whittetll2002h . An alterna- 
tive explanation is that the band strengths change due to interac- 
tion of H2O molecules with other constituents in the ice. In both 
high-mass and low-mass star forming regions, CO is - together 
with CO2 - the most dominant species that could mix w ith H2O. 
In a recent study on H20:C02 ices, IQberg et al ] (l2007h showed 
indeed significant band strength differences between pure and 
mixed H2O ices. The present study extends this work to CO con- 
taining water ice. 

C O accretes onto dust g rains around 20 K (San dford et al.l 

1988; Acharyya et al.l 12007) and plays a key role in solid state 
astrochemical processes, e.g., as a starting point in hydro- 
genation reactions t hat result in the formation of formalde- 
hyde and methanol ( Watanabe & Kouchil 20021; iHiraoka et al.l 
I2002t IWatanabe et al l 12004 iFuchs et al.l 120071 A strong ab- 
so rption centered arou nd 2139 cm was assigned to solid CO 
by ISoifer et al] dl979l) . again following thorough laboratory in- 
frared work. Further efforts in the laboratory have shown that 
CO molecules can be intimately mixed, either with molecules 
that possess the ability to form hydrogen bonds, such as H2O, 
NH3 and CH3OH - often referred to as "polar" ices - or with 
molecules that can only participate in a van der Waals type of 
bond, such as CO itself, CO2 and possibly N2 and O2 - so- 
called "non-polar" ices. In laboratory mixtures with H2O and 
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CO, the two forms are distinguished spectroscopically; the dou- 
ble Gaussian peak structure for the CO stretch fundamental can 
be decomposed in Gaussian profiles at 4.647 fim (2152 cirT 1 ) 
and 4.675 /im (2139 cm" 1 ), attributed to the polar and non-polar 
comp onent, respectively (ISandford et al.lll988t iJenniskens et alj 
1 1995b . On the contrary, pure CO measured in the laboratory ex- 
hibits a single Lorentzian band, which is located around 2139 
cm" 1 . This Lorentzian absorption profile can be further decom- 
posed into three Lorentzi an components cen tered around 2138.7, 
2139.7 and 2141.5 cm 1 ( Fraser et alj|2007h . 

In astronomical spectra, the 2139 cm" 1 feature has been con- 
sidered as an indicator of CO in H2O poor ic e, and the 2136 cm" 1 
feature as CO in H2O rich environments (Tielens e t"al]|199"Tl) . 
More recently it was found that the astronomical CO profiles can 
be decomposed into three components at 2136.5 cm" 1 , 2139.9 
cm" 1 and 2143.7 cm" 1 , with the 2139.9 cm" 1 feature ascribed 
to pure CO ice, and the 2143.7 cm" 1 feature ascribed to the lon- 
gitudinal optical (LO) component of the vibrat i onal transition 
in pur e crystalline CO (Pontoppid an et al.l l2003). Booger fet al.l 
(2002) proposed that the astronomically observed peak at 2143 
cm" 1 can originate from CO:CQ2 mixtures, but this i dentifica- 
tion is still controversial ( Ivan Broekhuizen et alj|2006l) . The as- 
signment of the 2136.5 cm" 1 feature in these phenomenologi- 
cal fits remains unclear. It should be noted that laboratory and 
astronomical data differ slightly in peak position, largely due 
to the fact that grain shape effects play a role for abundant ice 
molecules like CO and H2O. 

Recently, elaborate laboratory work and ab initio calcula- 
tions on mixtures of CO and H2O have shown that the absorp- 
tion around 2152 cm" 1 res ults from CO being bo und to the dan- 
gling OH site in H 2 ice dAl-Halabi et alj|2004 . Surprisingly 
enough, this absorpt ion has never been obse rved in the interstel- 
lar medium (e.g. iPonto ppidan et al. 2003). The non-detection 
of this feature has been explained by other molecules block- 
ing the dangling OH site, which is therefore unavailable to CO. 
An extension of this explanation is that the binding sites are 
originally populated by CO, but that this ha s been processed t o 
other molecules, such as CO2 or methanol (Fraser et al. 2004). 
Furthermore, it has been shown that the number of dangling OH 
sites decreases upon ion irradiation, which in turn results in a 
reduction of the integrated intensity of the 2152 cm" 1 feature 
dPalumboll2006l and references therein). The 2136-2139 cm" 1 
feature is a scribed to CO bound t o fully hydrogen bonded water 
molecules (lAl-Halabi et al.ll2004l) . 

Since CO and H2O are among the most abundant molecules 
in the interstellar medium, mixed CO and H2O ices have 
been subject to many experimental and theor e tical studies 
(e.g. iJiang et all 119751; iHagen & Tielens! Il98lt lHagen et al l 
1983; Al-Hala bi et al.l 120041; iFraser et alj 120051 For example, 
the behavior of the 2136-2139 cm" 1 CO stretching band has 
been quantitatively studied as a function of temperature and 
its band width and position have been studied as a function 
of H2O concentrati on in binary mixtures , but containing only 
up to 25% of CO dSchmitt et al.l Il989allbh . Furthermore, wa- 
ter clusters have been studied in a matrix of CO molecules 
with a ratio of 1:200 H20:CO. This has resulted in a ten- 
tative assignment of H2O monomers and dimers and the 
conc lusion that H2O forms a bifurcated dimer structure in 
CO (Hagen & Tielens! [T98 1 ). Other studies have focussed on 
Temperature Programmed Desorption (TPD) combined with 
Reflection Absorption Infrared Spectroscopy (RAIRS) of mixed 
and layered CO/H2O systems, enhancing great ly our knowl- 
edge o n their structures and phase transitions (Collings et al. 
2003a b). Nevertheless, a full quantitative and systematic study 



on the behavior of H2O in CO ice, and vice versa, with straight 
applications to astronomical spectra, is lacking in the literature. 
This is the topic of the present work. 

The desorption temperatures of CO and H2O differ by as 
much as 145 K under laboratory conditions. However, H2O/CO 
ices are expected to play a role in astronomical environments 
at temperatures not only well belo w the desorption temperature 
of CO at 20 K dFuchs et alj|2007l) . but also well above the des- 
orption temperature of p ure CO ice, since CO can be trapped 
in the pores of H2O ice (Collin gs et alj|2003al) . Thus far, both 
species have been observed together in lines of sight. It is often 
concluded from the non-detection of the 2152 cm" 1 feature that 
H2O and CO are not intimately mixed in interstellar ices. On the 
other hand, in some lines of sight CO is trapped in pores of a host 
matri x, as evidenced by the d etection of the 2136 cm" 1 CO fea- 
ture (Pontoppi dan et al1l2003l) . It is plausible that this trapping 
results from heating of a mixture of CO and a host molecule. 
Accordingly, we have also performed some experiments as a 
function of temperature. 

In this work, the effect of CO on the H2O vibrational fun- 
damentals is comp ared to the effect of CO2 on these modes, as 
studied recently by Obe rg et ail J2007). A comparison between 
the H2O bending mode characteristics in CO and CO2 contain- 
ing ices illustrates the sensitivity of this mode to the molecular 
environment. In addition, this work provides a unique laboratory 
tool for investigating the amount of CO mixed with water. 

The outline of this manuscript is as follows. In Sect. 2 the ex- 
perimental setup is described and the data analysis is explained. 
Section 3 is dedicated to the influence of CO on the water vibra- 
tional modes, as well as the influence of water on the CO bands. 
In Sect. 4, the astrophysical relevance is discussed and the con- 
clusions are summarized in Sect. 5. 

2. Experiment and data analysis 

The experimental se tup used for the measu rements has been de- 
scribed in detail in Gerakine s et alj d 1995b . It consists of a high 
vacuum setup (« 10" 7 Torr) in which ices are grown on a Csl 
window at a temperature of 15 K. The window is cooled down 
by a closed cycle He refrigerator and the sample temperature 
is controlled by resistive heating. A Fourier Transform InfraRed 
(FTIR) spectrometer is used to record ice spectra in transmission 
mode from 4000 to 400 cm" 1 (2.5-25 //m) with a resolution of 1 
cm" 1 . 

The sample gas mixtures are prepared in glass bulbs, using 
a glass vacuum manifold. The bulbs are filled to a total pressure 
of 10 mbar, which is well below the water vapor pressure. The 
base pressure of the manifold is better than 10" 4 mbar, result- 
ing in negligible contamination levels. A sample of CO (Praxair 
99,999%) is used without further purification. Deionized water, 
further purified by three freeze pump thaw cycles, is used for 
the H20:CO mixtures. Mixtures with different ratios H20:CO 
are prepared in the vacuum manifold and the resulting deposi- 
tions are listed in Table Q] The growth rate onto the ice is de- 
termined by setting the exposure to ~10 16 molecules cm" 2 s" 1 . 
Assuming a monolayer surface coverage of 10 15 molecules-cm" 2 
and a sticking probability of 1, this results in a growth rate of 10 
L-s" 1 (L = Langmuir). In the experiments where the effect of 
CO on the water ice vibrational modes is investigated, the water 
exposure has been kept constant with about 3000 L of water ice 
for the different mixtures to facilitate a one-on-one comparison 
between all samples. In the experiments where the effect of the 
H2O on the CO modes has been investigated, the total amount of 
deposited CO is kept constant at 3000 L (Table[TJ. 
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Table 1. Ice mixtures and resulting deposition thicknesses used 
in this work. Column A denotes the molecule of which the de- 
posited amount is kept constant, and column B indicates the 
molecule that is mixed in. The first series is used for determin- 
ing the effect of CO on the H2O band strengths and profiles. The 
second series is used to determine the effects of H2O on the CO 
band strengths and profiles. 
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10000 
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" 1 L (Langmuir) = 1 x 10 6 Torr s ~ 1 monolayer of molecules 

Three independent measurements are performed for 1:1 
H20:CO mixtures. These measurements allow for an estimate of 
the error in the experiment due to mixing of the gas, deposition 
of the sample and other errors that may occur. A conservative 
error of ^10% on the mixing ratios is deduced from these ex- 
periments. Additionally, two test measurements are performed 
for samples of 1000 and 10000 L, to check for layer thickness 
dependencies (TableQ}. 

The infrared transmission spectra are taken using a Biorad 
FTS40 spectrometer. A total of 256 spectra are acquired and 
averaged for each sample measurement. Background measure- 
ments are performed prior to each of the measurements and are 
subtracted to reduce the noise level. The spectra are further pro- 
cessed using IDL (Interactive Data Language) in order to flatten 
the baseline. This is done by fitting a second order polynomial 
through a set of five points, which are visually chosen well away 
from absorption features. The data reduction does not lead to a 
distortion of the absorption profiles. 

The absorption band strengths for the thre e modes of pure 
water ice at 15 K are well known from literature (Gerakine s et al.l 
11995b . The adopted values are 2xl0~ 16 , 1.2xl0~ 17 and 
3.1xl0~ 17 cm molecule -1 for the stretching (v s tretch=3279 cm -1 
or 3.05 /mi), bending (vb en d=1655 crrT 1 or 6.04 fim) and libra- 
tion mode (vi;b=780 cirT 1 or 12.8 //m), respectively (see FigJTJ. 
Calculating the integrated absorption bands over the intervals 
listed in Table [2] for the mixtures and scaling them to the inte- 
grated band strength for pure water ice, allows for a deduction 
of the band strengths for the water ice bands in the mixture via: 

X^band 
/H 2 0:CO=l:xX- (1) 

lnd Jband^ 

where AS™.., is the calculated band strength for the vibra- 
tional water mode in the l:x mixture, 7h,0:CO=1:x its integrated 
area, A^ al ^ the band strengths for the water modes as available 

from literature and f 7h,o the integrated area under the vibra- 
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Fig. 1. A typical baseline and background corrected ice spectrum 
for a H20:CO = 1:1 mixture. The measurement is performed at 
15 K. The vibrational modes in the H20:CO ice are indicated. 



Table 2. The measured peak positions and the integration bounds 
in cirT 1 used to compute the integrated intensities of the H2O 
bands. The values between brackets indicate the /mi values. 









Integration bounds 


Species 


Assignment 


Peak 


Lower 


Upper 


H 2 


^libration 


780 (12.8) 


500 (20.0) 


1100 (9.09) 




y bend 


1655 (6.04) 


1100 (9.09) 


1900 (5.26) 




y stretch 


3279 (3.05) 


3000 (3.33) 


3600 (2.78) 




Vfree OH 


3655 (2.73) 


3600 (2.78) 


3730 (2.68) 


ccr 


V' stretch' 


2139 (4.68) 


2120 (4.72) 


2170 (4.61) 




V' polar' 


2152 (4.65) 


2120 (4.72) 


2170 (4.61) 



" The integrations for the two CO bands are performed using 
Gaussian fits. 



tional mode for the pure water sample. The free OH stretching 
mode, 'the fourth band', is scaled to the stretching mode for pure 
water since this absorption is absent in the spectrum of pure H2O 
ice. 

Integration limits used throughout the experiment are listed 
in Table [2] Integrated areas relative to the integrated area of 
the pure water stretching mode, A/ A pure H 2 stretch, are in- 
vestigated as a function of CO concentration. For the sam- 
ple with the most mixed in CO, i.e., the 1:4 H20:CO mix- 
ture, an analysis in terms of cluster formation is given. In 
addition, the influence of temperature on the water stretch- 
ing mode is studied. The measured spectra for the H20:CO 
mixtures are available online at the Leiden ice database: 
http://www.strw.leidenuniv.nl/~lab/databases 

3. Results 

3. 1. Influence of CO on water bands 

In Fig. [2] the four H2O ice fu ndamentals are sho wn for different 
compositions. Similar to CO2 (Oberg et al. 2007), CO has a clear 
influence on the water ice absorption bands compared to the pure 
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Fig. 2. Combined spectra of the four modes for water ice for the six measured compositions (see Table [T]), ranging from pure water 
ice (bottom) to a 1:4 H20:CO mixture (top figures). The spectra are taken at a temperature of 15 K. Note that the wavelength ranges 
for separate modes are different. The small structures on the libration mode are experimental artifacts. 
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Fig. 3. The integrated intensity of the water vibration modes rel- 
ative to the integrated intensity of the pure water stretch mode 
plotted versus the concentration of CO in the sample ice. It 
should be noted that the plots for the bending, free OH and libra- 
tion mode have been multiplied by a factor of ten to facilitate the 
display and that the stretch mode is the one that is actually most 
dependent on the concentration. The four water modes show to 
first order a linear dependence on the CO concentration. The es- 
timated error in the data amounts to 10%. 



H2O ice. This effect is different for each of the four bands. The 
bulk stretch mode is most strongly affected; the band strength for 
this mode decreases by more than a factor of 2 when the CO frac- 
tion is raised from 0% to 80%. The band strength of the free OH 
stretch, which is absent when no CO is mixed in, is greatly en- 
hanced with concentration. The libration mode gradually looses 
intensity when the amount of CO in the ice is increased and the 
peak of the absorption band shifts to lower energy. The inte- 
grated areas of the four water modes are scaled to the pure wa- 
ter stretching mode and plotted versus the CO concentration in 
Fig. [3] A linear function A e ft = a- [CO] + b is fitted through the 
data points of the four water modes. The a coefficient indicates 
the strength of increase/decrease of the band strength, and the 
b coefficient indicates the band strength of water relative to the 
pure stretching mode when no impurities are mixed in. There ex- 
ists some deviation between the data points and the fit function, 
which is most probably due to the deposition accuracy, but this 
deviation is within the experimental error of 10%. A clear trend 
in all four modes is observed. In Table[3] the linear fit coefficients 
are listed for the H20:CO binary mixtures. The linear coefficient 
for the H2O stretching mode is highest and negative, indicating 
the strongest decrease in band strength. The free OH stretching 
mode has a positive linear coefficient indicating that this is the 
only mode to increase in intensity upon CO increase. A com- 
parison with recently obtained data for H20:C02 ices shows the 
same trend. Apart from the bending mode, all effects are more 
pronounced in the CO2 mixtures, i.e., the absolute values of the 
a coefficient are larger by a factor of 1.3 - 2. This is related to 
the actual interactions in the ice and work is in progress to study 
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such effects in more detail for a large number of species from a 
chemical physics perspective. 

The free OH mode, the water stretching mode and the water 
bending mode start showing substructure superposed onto the 
bulk absorption profile upon increase of the fraction of CO in 
the ice mixture (Fig. @}. The absorptions of the bulk modes are 
still clearly apparent beneath the substructure. For the stretch- 
ing mode this absorption shifts from 3279 crrT 1 to a higher 
wavenumber of 3300 crrT 1 . The peak absorption for the bulk 
bending mode shifts from 1655 cirT 1 for pure FbO to 1635 cm 1 
for the 1:4 H20:CO mixture. The libration mode is also red- 
shifted upon CO concentration increase. For the pure H2O ice 
spectrum this mode is located at 780 cm -1 , while for the high- 
est partner concentration it appears at 705 cm"'. The free OH 
stretch gradually increases in frequency. The peak absorption 
shifts from 3636 cm -1 for the 1:0.25 to 3655 cm" 1 for the 1:4 
H20:CO mixture, which corresponds to a blue shift of 19 cm"'. 

The substructure, which is superimposed on the bulk stretch- 
ing mode around 3300 cm" 1 , has been previously assigned to 
(H2O),, water clusters in the ice. These assignments are based 
on matrix spectroscopic d ata of H2O in a matrix of N2 + O2 
(75:25) dOhno et al.ll2.Q05b . Comparable data for H2O in a ma- 
trix of CO are not available. Nevertheless, in Fig. |4]it is shown 
that an excellent fit is obtained when peak positions from matrix 
spectroscopic data are used, assuming that the peak positions are 
shifted, which is indicative for the difference in interaction be- 
tween H2O and CO compared to H2O and N2 + O2. For each of 
the contributions, the peak position, bandwidths (Full-Width-at- 
Half-Maximum, FWHM) and integrated area, are sum marized in 
Table El and compared to previous results obtained bv lOhno et alj 
(120051 Some absorptions are red-shifted and other absorptions 
are blue-shifted compared to absorptions of H2O clusters in a 
N2 + O2 matrix. Note that the relative H2O concentrations in 
the present work are sub s tantia lly higher than in the matrix ex- 
periments bv lOhno et ail 12005). As a consequence, larger H2O 
clusters are more pronounced in our spectra. 

Additionally, the temperature of the sample plays a role on 
the band strengths. A clear effect is encountered when the tem- 
perature is slightly increased, i.e., from 15 to 25 K, close to 
the desorption temperature of CO. At this temperature, the CO 
molecules in the matrix start gaining enough energy to become 
mobile. The mobility in the matrix allows water clusters to find 
partners for hydrogen bindings and to reorganize themselves to 
form a stronger bulk hydrogen bonded network, as indicated by 
the increased bulk stretch mode band strength and decreasing in- 
tensity of the substructure. Figure [5] nicely shows the transition 
from H2O clusters embedded in a matrix of CO to the formation 
of a bulk hydrogen bonded network. 

Figure [6] shows the different effect of CO and CO2 on the 
H 2 bending mode in ice for the 1 :4 H 2 0:C0 2 and 1 :4 H 2 0:CO 
mixtures. The intensity ratio of the main peaks is actually re- 
versed in the two mixtures. The differences between the two 
mixtures start showing up from a mixing ratio of 1:0.5 H20:X 
and become more pronounced for higher CO and CO2 concen- 
trations. In addition, the CO mixtures exhibit a stronger broad 
underlying feature, which is visualized by the Gaussian fit in 
Fig. [6] In other words, a detailed study of the H2O bending 
and stretching modes may provide additional information on 
whether CO or CO2 dominates in the ice. The free OH stretch- 
ing mode is also affected differently by the two molecules. In 
the H20:C02 mixtures, the mode is more shifted to higher 
wavenumbers. For the 1 : 1 H20:C02 mixture the peak position is 
3661 cm" 1 , compared to 3635 cm" 1 for the H20:CO mixtures. 
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Fig. 4. Measured (indicated in black) and fitted spectrum (indi- 
cated in red) of a 1:4 H20:CO mixture in the frequency range of 
the H2O bulk stretch and free OH stretch modes at a temperature 
of 15K. The water stretching mode clearly shows substructure. 
An excellent fit is obtained by a superposition of the bulk stretch 
mode and Gaussian functions representing s maller H2O cluster s 
(indicated in green) in the matrix material (iOhno et aTJl2005h. 
The fi tted peak positions and the positions from Ohn o et all 
(2005) are listed in Table [4] Please refer to the online version 
for the color coding. 



Thicker and thinner layers of the mixtures have been mea- 
sured to check for thickness dependence. We conclude that 
within our experimental error limit, ice thickness does not play a 
significant role in the behavior of the relative band strengths. 
This conclusion is supported by the observation that identi- 
cal mixing ratios in the two measurement series (H2O/CO and 
CO/H2O) show the same (scaled) spectroscopic behavior for dif- 
ferent total ice thicknesses (Table[T]i. 

3.2. Influence on the CO band 

Section 3.1 shows that mixing CO into a water mixture affects 
the band strengths of the water vibrational modes. Reciprocally, 
the CO stretch mode is also altered when water is added to CO 
ice. This is seen in the experiments where the amount of de- 
posited CO is kept constant (TableQ]). When water is mixed into 
the CO ice, the absorption changes from a Lorentzian profile 
to a Gaussian profile. Furthermore, the second CO absorption 
at 21 52 cm" 1 , ascribed to CO bound to the H2O dangling OH 
sites dFraser et alj 12004). manifests itself as a Gaussian profile 
and increases in integrated intensity upon increase of water con- 
centration. The transition from the pure CO Lorentzian shaped 
profile to two Gaussian shaped profiles for the mixed CO:H20 
ices is illustrated in Fig. [7] We have not further decomposed 
the 2139 cm"' com ponent. More information is available from 
lFraseretal.1 (120071) . 

Gaussian fits for both CO absorption components are made 
for the range of mixtures as listed in Table Q] at a temperature of 
15 K. One typical fit is shown in Fig. [8] With this example, it is 
demonstrated that excellent fits are obtained by using the fit pa- 
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Table 3. Resulting linear fit coefficients for the H20:CO mixtures. The coefficients indicate the strength of the interaction between 
CO and the H2O host molecules in the matrix for mixtures that are deposited at a temperature of 15 K. The corresponding values 
for H20:CC>2 ice mixtures are listed for a comparison. 







Linear Coefficients 




mixture 


H 2 mode 


constant (b ) 
[1CT 16 cm molecule -1 ] 


slope (a ) 
[1(T 19 cm molecule -1 ] 


R 2 


H 2 0:CO 


^libration 


0.30±0.02 


-2.1±0.4 


0.93 




Vbcnd 


0.13±0.02 


-1.0±0.3 


0.84 




y stretch 


2.0±0.1 


-16±3 


0.95 




v fiee OH 





1.2+0.1 


0.99 


H 2 0:C0 2 a 


Vlibration 


0.32±0.02 


-3.2±0.4 


0.99 




v bcnd 


0.14±0.01 


-0.5±0.2 


0.81 




"^stretch 


2.1+0.1 


-22±2 


0.99 




Vficc OH 





1.62±0.07 


0.99 



" From Obers et al. (2007 



Table 4. Line positions, FWHMs and integrated areas of the Gaussian functions fitted to the 1:4 H20:CO water bulk stretch and 
free OH stretch spectru m. The assignment is based on both Density Functional Theory (DFT) calculations and experimental values 
from lOhno etail (E005) obtained in a N2/O2 matrix. 



Mode 


Position" 


Position'' 


A Position 6 


FWHM" 


Area'' 




(cm -1 ) 


(cm" 1 ) 


(cm- 1 ) 


(cm- 1 ) 


(a.u. • cm 1 ) 


Free OH (ring) cyclic-pentamer, cyclic-trimer free OH, clusters hexamer 


3688 


3674 


15 


15 


0.371 


Asymmetric stretch H-acceptor dimer 


3715 


3693 


22 


7 


0.040 


Asymmetric stretch monomer 


3715 


3706 


9 


5 


0.043 


Symmetric stretch monomer 


3635 


3608 


27 


5 


0.074 


Symmetric stretch H-acceptor dimer 


3629 


3596 


33 


15 


0.078 


Bulk free OH stretch 




3655 




19 


1.57 


Bulk free OH stretch 




3633 




25 


0.609 


H-bonded OH stretch "chair" hexamer 


3330 


3331 


-1 


43 


0.987 


H-bonded OH stretch "cage" hexamer 


3224 


3226 


-2 


34 


0.626 


H-bonded OH stretch "prism" hexamer 


3161 


3175 


-14 


40 


0.119 


H-bonded OH stretch pentamer 


3368 


3364 


4 


29 


0.491 


H-bonded OH stretch "bookl, cage" hexamer 


3450 


3420 


30 


86 


2.62 


Bulk stretching mode 




3300 




250 


11.5 


H-bonded OH stretch ring cyclic trimer 


3507 


3506 


1 


42 


0.894 



" Ohno et al. (2005), " This work. 



rameters as listed in Table [5] The behavior of the integrated area 
of the 2152 cirr 1 component compared with the total integrated 
CO absorption is plotted as a function of CO and H2O concentra- 
tion in Fig. [TO] The integrated area for the 2152 cirr 1 component 
decreases with increasing CO content, and it is undetectable for 
a pure CO ice. A second order polynomial describes how the 
polar component behaves with respect to the CO concentration 
[x] or water concentration [100 - x] in the ice over the interval 
spanning from 20% CO up to a pure CO ice. Note that the total 
amount of deposited CO is kept constant. The coefficients of the 
second order polynomial of the form y = a-x 2 + b- x + c are 
a=-0.005, b-0.23 and c=26.6. For a decreasing amount of water 
in the sample, the peak position of the 2152 cm-' absorption fea- 
ture is most strongly affected and decreases gradually to lower 
wavenumbers, until it reaches 2148 cirr 1 for the 1:0.25 CttHbO 
mixture. The FWHM of this band is also affected. It starts at a 
width of 10.5 cm 1 for the 1:4 CO:H20 mixture and decreases 
to 7.5 cm 1 for the 1:0.25 CO:H20 mixture. The position of the 
main absorption feature at 2139 cm -1 is only slightly affected 
by increasing the amount of water in the sample and decreases 
by 1.3 cm 1 when going from pure CO ice to the 1:4 CO:H20 
mixture, i.e., a shift towards the 2136 cirr 1 feature is observed 
upon dilution (see Fig. |9). Its position is expected to shift even 
more, to 2136 cm -1 , for H2O concentrations above 80%. The 
FWHM of this absorption feature decreases from 8 cirr 1 for 
the 1:4 CO:H 2 to 5 cm 1 for the 1:0.25 CO:H 2 sample. The 



Lorentzian peak profile of the pure CO absorption exhibits an 
even smaller FWHM of 2 cm-'. An overview of the changes in 
peak position, FWHM and integrated intensity is given in Table 


The position of the 2139 cm 1 absorption feature is also 
strongly dependent on the temperature of the ice, as illustrated 
in Fig. [9] The majority of the CO will desorb as the ice is heated 
above the CO desorption temperature. The remaining CO shows 
an absorption that is shifted toward 2135 cm 1 . Thus, a shift 
from 2139 to 2136 cirr 1 occurs both by mixing with significant 
amounts of H2O and by heating of CO:H20 above 40 K, even for 
mixtures with modest amounts of H2O. From Fig. [9] it becomes 
clear that the latter effect (i.e. heating) is the more critical one. 
One should note that the laboratory data presented here cannot 
be compared one-to-one with the observational data because of 
shifts caused by grain shape effects. 

4. Discussion 

CO and H2O ice abundances are often derived using the well 
known constants from the literature. The experimental work pre- 
sented here shows that band strengths deduced from pure ices 
cannot be used to derive column densities in interstellar ices 
without further knowledge on the environmental conditions in 
the ice. Concentrations of CO and CO2 ice as high as 15 and 
21%, respectively, relative to H2O, have been reported towards 



Jordy Bouwman et al.: Band profiles and band strengths in mixed H20:CO ices 



7 



2140.5 -■ 

2140.0- 

2139.5- 

'£ 2139.0- 
o 

2138.5- 

O 

£ 2138.0- 

(0 

O 

a- 2137.5- 
J£ 

CO 

o 2137.0- 
Q. 

2136.5- 
2136.0- 
2135.5 4- 



- r- 

20 



- r- 

40 



- 1 — 

60 



— I — 

80 



2140.5 -■ 
2140.0- 
2139.5- 
2139.0- 
2138.5- 
2138.0- 
2137.5- 
2137.0- 
2136.5- 
2136.0- 
2135.5- - 



60 



80 



Water concentration / % 



Temperature / K 



Fig. 9. Change in peak position of the 2139 cm 1 CO absorption feature as a function of both H2O concentration and temperature 
for the 1:4 H20:CO mixture. The dotted lines indicate the band positions measured for CO in the interstellar medium. 



Wavelength / urn 
2.8 3 3.2 

0.12 - — , 1 , 1 , ,- 

30 K 




Fig. 5. Temperature dependence of larger clusters of water 
molecules in a matrix of H20:CO = 1:4. For increasing temper- 
ature the substructure gives way to the bulk stretch mode when 
CO evaporates. 

GL7009S dKeane et al.ll200Th . If we assume a polar fraction of 
75%, as reported towards W33A, and that the polar fraction of 
both CO and CO2 ice are in close contact with the water, the 
band strength for the H2O bending mode is reduced by a fac- 
tor of ~ 1.25. In other words, the band strength for H2O will be 
smaller and hence the column density of H2O will be underesti- 
mated if the laboratory data from pure H2O ice are adopted from 
the literature. 

Vice versa, for deriving the column density of CO ice one 
should also consider the influence of other molecules in its vicin- 



Wavelength / jim 
5.8 5.9 6 6.1 6.2 6.3 6.4 




1 1 1 1 1 1 1 

1750 1700 1650 1600 1550 



Wavenumber / cm" 

Fig. 6. The spectral differences in the wate r bending mode pro- 
file for a 1:4 H 2 0:C0 2 ice mixture (top) dOberg et alj2007l) and 
a 1:4 H20:CO mixture (bottom). The black bold (overall) spec- 
tra indicate the laboratory spectra and coincide with the fitted 
spectra (red dotted) consisting of the three Gaussian curves, in- 
dicated by the blue dotted lines. (See online version) 



ity. The integrated absorbance of CO, however, is less strongly 
influenced by the presence of H2O. The CO absorption decreases 
in absolute intensity when water is mixed in, but the total in- 
tegrated band strength is compensated by a broadening of the 
absorption and the appearance of the 2152 crrr 1 component, at 
least in the laboratory spectra. Upon heating, the CO molecules 
become trapped into pores, as indicated by the 2136 crrr 1 fea- 
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Table 5. Lorentzian and Gaussian fit parameters for the CO 
stretching mode for a constant amount of CO in ice mixtures 
ranging from 100% CO to a 1:4 CO:H 2 mixture. 



Composition 


Position 


FWHM 


Area 




(cm- 1 ) 


(cm- 1 ) 


(a.u. • cm" 1 ) 


Pure CCT 


2138.8 


2.2 


1.9 


1:0.25 CO:H 2 


2138.2 


5.0 


1.31 




2147.5 


7.5 


0.23 


1:0.5 CO:H 2 


2138.0 


5.8 


1.40 




2148.1 


7.8 


0.34 


1:1 CO:H 2 


2138.2 


6.5 


1.28 




2148.3 


9.0 


0.44 


l:2CO:H 2 


2137.9 


7.6 


1.24 




2149.5 


9.8 


0.49 


l:4CO:H 2 


2137.5 


7.9 


1.38 




2150 


10.5 


0.58 



" Lorentzian profile 



Fig. 7. Illustration of the behavior of the CO stretch fundamen- 
tal upon increase of the concentration of H 2 in CO. The total 
integrated band strength remains unchanged within the experi- 
mental error, although the maximum intensity of the absorption 
decreases strongly. The y-axis is cut off for the pure CO mode to 
make a clearer distinction between 'non-polar' and 'polar' com- 
ponents of the CO absorption for the CO:H20 mixtures. 



4.62 



Wavelength / |xm 
4.64 4.66 4.68 



4.70 



0.15- 




2170 



2160 



2150 2140 
Wavenumber / cm" 1 



2130 



2120 



Fig. 8. Gaussian fit of the CO stretch mode in a CO:H 2 = 
1:0.25 mixture. The black bold spectrum is the measured lab- 
oratory spectrum that is reproduced (red line) by adding the two 
Gaussian components centered around 2138.2 cirr 1 and 2147.5 
cirr 1 (green dotted lines). These are attributed to CO in a 'non- 
polar' and 'polar' environment. (See online version) 



ture. The abso rption strength of this band is very sensitive to 
temperature (Schmitt et al. 1989a). The CO absorption profile is 
indicative of the amount of water that is mixed. The percentage 
of CO mixed into the H 2 ice is derived in the laboratory from 
the ratio between the 2152 cirr 1 feature and the total CO absorp- 




Fig. 10. Absorbance of the 2152 cm 1 component relative to the 
total CO absorption as a function of CO concentration. 



tion (Fig. [TOt . This allows the derivation of the effective band 
strength A^ nd for H 2 in a mixture with CO using the linear 
model proposed in Sect. 3.1, and t hus an estimate of th e column 
density for H 2 via the equation (Sand ford et al.l fl988): 



J Tydv 

A eff 



(2) 



Temperature also affects the band strength as shown in Fig. 
[5] The applications of the model presented in Fig. [10] are re- 
stricted to deduction of column densities for ices with a temper- 
ature of about 15 K, which have not been thermal ly processed. It 
should also be noted, as recently demonstrated by Biss chop et al.l 
(2007), that in addition to binary mixtures, tertiary mixtures also 
have to be taken into account to compare laboratory data with 
astronomical spectra. 
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5. Conclusions 

Based on the experiments described in this manuscript and in 
recent work on H^OiCCh mixed ices, we draw the following 
conclusions regarding the interaction between H2O and CO in 
a solid environment: 

1. The general trend on the band strengths of the four vi- 
brational modes in water ices is similar for H^OtCO and 
H 2 0:C02 mixtures with increasing CO or CO2 concentra- 
tion. However, quantitative differences exist, reflecting dif- 
ferences in the strength of the interaction, which allow us to 
distinguish between CO and CO2 in H2O ice, explicitly as- 
suming that the main constituents of the ice are H2O, CO2 
and CO. 

2. The position of the water free OH stretching mode is partic- 
ularly indicative of the molecule that is interacting in the ma- 
trix (CO vs CO2), again under the assumption that we only 
consider H2O/CO2/CO ices. The peak position of this mode 
is 26 cm -1 red-shifted for a 1:1 H20:CO ice mixture com- 
pared to a 1:1 H20:C02 ice mixture. 

3. In addition, the water bending mode is indicative of the 
molecule, i.e., CO or CO2, that is interacting with the water 
ice. The relative integrated intensity of the Gaussian compo- 
nents reveals whether CO2 or CO is mixed into the H2O ice. 
The same restriction mentioned in conclusions 1 and 2 also 
applies here. 

4. Upon increasing the relative amount of CO in the mixture, 
a clear substructure begins to reveal itself in the bending, 
free OH stretch and bulk stretching mode. The arising sub- 
structure indicates the onset of H2O cluster formation in the 
H20:CO ice. An assignment of the clusters was possible fol- 
lowing matrix isolation spectroscopy. 

5. The substructure on the stretching mode quickly gives way 
to the bulk water mode when the temperature is increased 
close to the desorption temperature of CO. This can easily 
be depicted by CO molecules becoming mobile and hence 
allowing single water molecules and larger water clusters to 
find partners for bulk hydrogen bonding. 

6. The ratio 2152 cm -1 / total integrated CO absorption inten- 
sity is a tracer of the amount of CO that is mixed into the 
laboratory water ice, or vice versa. In astronomical spectra 
this band has not been observed. 

7. H2O column densities derived from astronomical spectra can 
easily be underestimated by as much as 25% when environ- 
mental influences, i.e., CO or CO2 presence, are not taken 
into account. 

The present systematic study of CO:H20 ice, together with 
recent work on C02:H20 ice, provide the tools with which to 
estimate the mixing ratios of the three most abundant molecules 
in interstellar ices. 
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